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Abstract 

The exact ground state of the many-body Schrodinger equation for N bosons on a one- 
dimensional ring interacting via pairwise (^-function interaction is presented for up to fifty particles. 
The solutions are obtained by solving Lieb and Liniger's system of coupled transcendental equa- 
tions for finite N. The ground state energies for repulsive and attractive interaction are shown to 
be smoothly connected at the point of zero interaction strength, implying that the Bethe-ansatz 
can be used also for attractive interaction for all cases studied. For repulsive interaction the exact 
energies are compared to (i) Lieb and Liniger's thermodynamic limit solution and (ii) the Tonks- 
Girardeau gas limit. It is found that the energy of the thermodynamic limit solution can differ 
substantially from that of the exact solution for finite N when the interaction is weak or when iV 
is small. A simple relation between the Tonks- Girardeau gas limit and the solution for finite inter- 
action strength is revealed. For attractive interaction we find that the true ground state energy is 
given to a good approximation by the energy of the system of N attractive bosons on an infinite 
line, provided the interaction is stronger than the critical interaction strength of mean-field theory. 
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I. INTRODUCTION 



The recent experimental realization of quasi one-dimensional Bose-Einstein condensates 



QflfiQfl 



has made the theoretical description of these systems a very active area of 



research. In this context, the underlying equation is the many-body Schrodinger equation for 
N particles subject to two-body ^-function interaction and possibly an external potential 
^]. In general, solving this equation exactly is very difficult, and it is almost inevitable 
to introduce approximations. The most commonly used approximation for Bose-Einstein 
condensates is the Gross-Pitaevskii (GP) approximation 0,^^. The GP approximation is 
a mean-field approximation that results in a comparatively simple nonlinear equation. This 
equation can be solved numerically and explicit analytical solutions have been found for some 
cases, see [ill 12, 3 and references therein. Despite its great success in the description of 
early experiments, the GP approximation suffers from various shortcomings. For instance, 
the solutions _of the GP eq uation may not possess the symmetry of the Hamiltonian of 

More recently, other approximations have been developed to 
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the problem 

Exactly solvable one-dimensional models are of interest by themselves and can be consid- 
ered as a research discipline of its own, see l^cj] for an overview. Moreover, they provide 
an invaluable testing ground for approximative methods. In fact, this is often the motivation 
for the research in this field. However, from an experimental point of view the assumption 
of one-dimensionality seems far-fetched at first sight. Therefore, it is very exciting to see 
the experimental realization of some of these model systems within reach. 

A very prominent exactly solvable model is Lieb and Liniger's system of N spinless 
ooint-like bosons on a one-dimensional ring interacting via pairwise 5-function interaction 



2l| . This is also the subject of the present work, and the proper definition of the problem is 



given in the next section. Lieb and Liniger's model is a generalization of Girardeau's gas of 
impenetrable bosons [3] to finite interaction strength. The impenetrable boson gas is also 
known as the Tonks-Girardeau (TG) gas, thereby including also the name of the inventor of 
the classical hard sphere gas 2^| . The TG gas is not less prominent than its finite interaction 
counterpart, and even though both models are now more than forty years old, they are still 

26. 



the subject of ongoing research, see 
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27, 



In their ground-breaking work Lieb and Liniger derived a system of iV — 1 coupled tran- 
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scendental equations that determine the exact iV-particle ground state of the problem. Lieb 
and Liniger solved this system explicitly for two particles, but then passed to the thermo- 
dynamic limit of the system which exists for repulsive interaction only. Surprisingly, in this 
limit the whole system of coupled transcendental equations can be approximated by a single 
Fredholm integral equation of the second kind. Lieb and Liniger solved this integral equa- 
tion already in their initial work. Moreover, they proved that it has an analytic solution 
for any interaction strength truly greater than zero. For weak interaction, Bogoliubov's 
perturbation theory agrees well with this solution and for strong interaction its energy con- 
verges to that of the TG gas in the thermodynamic limit 22j. In a subsequent paper Lieb 
also derived the excitation spectrum of the thermodynamic limit solution [29]. Lieb and 
Liniger's thermodynamic limit solution has also been used to describe interacting bosons 
in more general one-dimensional trapping potentials than just a one dimensional ring. For 
example, by assuming that the thermodynamic limit approximation is locally valid an d by 
employing a hydrodynamic approach, bosons in cigar shaped traps have been described [30] . 

However, the set of coupled transcendental equations which yield the exact solution of 
the finite N problem was not solved for N > 2 until 1998, neither for repulsive nor for 
attractive interaction. In 1998 Muga and Snider derived the whole spectrum of the three- 
particle problem j^| for attractive and repulsive interaction. Still, up to date there is no 
exact solution of the problem for N > 3. In the present work we would like to fill in this 
gap. We calculate the ground state for up to fifty particles. For repulsive interaction even 
for fifty particles the thermodynamic limit solution can deviate from the solution for finite 
N by as much as two percent, as we shall show. For attractive interaction we reveal a close 
relation to the system of N interacting bosons on an infinite line, see 14] and references 
therein. 

It has been proven that for repulsive interaction the wave functions of all states are of the 
Bethe-ansatz type This is not the case for attractive interactions Q,0]. However, 

Muga and Snider have shown that for three attractive bosons j^j] a complete set of states 
can probably be derived from a Bethe-ansatz. We will clarify this at a later stage. In the 
present work we give further support to the hypothesis that all states may be derived from 
a Bethe-ansatz for all particle numbers and all interaction strengths. 

As we have shown recently, the problem of iV bosons on a ring defies any accurate 
description by using direct diagonalization techniques, except for very weak interaction [16] . 
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In the present work, we present the exact ground state solution of the system of N bosons 
on a one-dimensional ring by solving Lieb and Liniger's system of coupled transcendental 
equations for up to fifty particles. We treat repulsive and attractive interactions alike and 
compare the results with other limiting cases of this system. Our approach can easily be 
extended to any particle number and also to excited states. 

The plan of the paper is as follows. In section |H] we define the problem, review the 
derivation of Lieb and Liniger's system of coupled transcendental equations and derive re- 
sults for weak attractive and repulsive interaction, by using first order perturbation theory. 
Section IHII addresses the problems related to the numerical solution of the set of coupled 
transcendental equations. The attractive case proved to be particularly delicate. We show 
that the Bethe-ansatz gives solutions for attractive and repulsive interaction that are con- 
tinuously connected to each other at the point of zero interaction for all particle numbers 
under consideration. In section HVl we present the exact energies for repulsive interaction for 
up to fifty particles. We compare these energies with the energies of first order perturbation 
theory, the TG gas and the thermodynamic limit. Furthermore, we present the differences 
of these limits to the finite N solution and reveal a surprisingly simple relation between the 
TG limit and our solution. We give explicit limits on the minimal number of particles and 
the size of the ring for the thermodynamic limit approximation to be reasonably accurate. 
In section El the results for attractive interaction are presented. The exact ground state en- 
ergies are found to converge to the energies of the one-dimensional problem of N attractive 
bosons on an infinite line, see Q] and references therein. We give a simple explanation for 
this behaviour in a mean-field picture and show that this energy is approached in the limit 
of very strong attractive interaction. To be concrete, we also calculate the minimal number 
of particles which is necessary such that the system on a ring of finite size can be approxi- 
mated by the system on an infinite line. Section IVII contains a summary of our results and 
a discussion of open questions. 
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II. THEORY 



A. The Schrodinger equation of the problem 

Our starting point is the stationary Schrodinger equation for N bosons in one dimension 
subject to two-body ^-function interaction and periodic boundary conditions: 

i=l y * i<j 

where the are the particle coordinates < yi < I and 

il)(y 1 ,...,y u ...,y N ) = ip(y 1 , ...,y t + I, ...,y N ), (2) 

for z = 1, JV. We divide equation ((TJ by and change from the dimensional to dimen- 
sionless coordinates Xi = j-L, where L is the new dimensionless length of the ring. The 
combined effect of these changes results in the equation 

- d 2 tb N 

i=l * i<j 

where now < X{ < L. The relation between the dimensionless and the dimensional 
quantities is given by the equations 

Vi = jl, (4) 

"srr m 

From (jSJ and (jSJ it follows that the relation between the dimensionless energies of two rings, 
one of length L and one of length L', is given by the equation 



Mc) = (7) 

L 2 



Attractive interactions are described by c < and repulsive by c > 0. 
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B. The Bethe-ansatz wave function 



The ^-function potential in (jSJ) is equivalent to a jump in the derivative of the wave 
function, wherever two particles touch [21| : 

~- 2ajj\ Xj=Xk . (8) 



fdlf) 


dip\ 


/dip 


dip\ 


\dXj 




\dx k 


^"^j J Xj=Xf c 



As long as c < oo the interaction potential allows the particles to go past each other, and 
any particle can be anywhere in coordinate space. However, since the particles are identical 
bosons the knowledge of the wave function in the "primary" region 

R P ■ < x x < x 2 < ... < x N < L, (9) 

namely 

ip(xi < x 2 < ... < x N ) (10) 

contains the full information. It has to be stressed that the wave function in R p does not 
have to be symmetric. Once the wave function in R p is known, the wave function in any 
other region corresponding to a different ordering of the coordinates is obtained simply 
by interchanging the particle labels in (fTU|) . This ensures the symmetry of the total wave 
function under particle exchange in the unrestricted coordinate space. In the region R p , 
equations Q and (JBJ) become 

N 8 2 

-£^> = £l(c#, (11) 



, dx] 
i=i 1 



for Xi 7^ Xj and 



d d 



dxj + i dxj 



^\x j+1 = Xj = ap\ Xj+1=Xj . (12) 

Since periodic boundary conditions are used, a displacement of L in any of the coordinates 
leaves the wave function unchanged. Applying a displacement of L to the wave function in 
R p yields 

^(0, x 2 , ...,x N ) = tp(x 2 , x 3 , x N -i, L) (13) 

and for the derivatives 

d d 

—ip{x,x 2 , ...,x N )\ x=0 = ^r?p(x 2 ,x 3 , ...,x N ,x)\ x=L . (14) 
ox ox 
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Two particles only interact when they are at the same point in space. This fact and (jlljl 
motivate the idea that the solution might be a product of plane waves. In fact, the Bethe- 
ansatz is just a more general version of this idea, namely a superposition of products of 
plane waves. It was first applied to spin chains but has been successfully used to solve 
numerous other one- dimensional problems exactly 20|. The Bethe-ansatz wave function for 



this problem is 



21] 



N 



ip(xi < x 2 < ... < a; at) = a {P) P Gxp(i kjXj) , (15) 

p j=l 

where the sum runs over all permutations P of the {xj}, and the a(P) are coefficients which 
are determined by the rule given below. It has been proven that the Bethe-ansatz gives 
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For the Bethe-ansatz to be 



all solutions of the problem for repulsive interaction 
valid, all fcj must be different from one another if c ^ 0, otherwise ip vanishes identically by 
means of (|12p. Only if c = 0, the equality of two kj does not imply a vanishing Bethe-ansatz 
wave function. However, there are certain critical c values at which the equality of two k~ 



does occur 



21 



[31 1 . implying that the form of the Bethe-ansatz wave function is no longer 
valid. This only reflects the fact that the normalization constant has not been included in 
the definition of the Bethe-ansatz wave function. The proper normalized wave function can 
be obtained by using the rule of l'Hospital As Muga and Snider have shown, for three 
particles all Bethe-ansatz solutions are continuously connected in /c-space at c = 0. It is 
therefore very likely that for three particles a complete set of states can be derived from the 
Bethe-ansatz in the sense that l'Hospital's rule is applied to obtain the wave function from 
the Bethe-ansatz at the critical c-values. As we shall show in section the Bethe-ansatz 
provides ground state solutions that are continuously connected in fc-space at c = 0, at least 
for as many as fifty particles. Therefore, we suspect that a complete set of states can be 
derived for all particle numbers from a Bethe-ansatz wave function in the sense described 
above. 



C. Lieb and Liniger's transcendental equations 

In the following we review the derivation of Lieb and Liniger's system of coupled tran- 
scendental equations as far as it is indispensable for our needs. However, we adopt Muga 
and Snider's approach and notation, since it allows to treat the repulsive and the attractive 
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case in a coherent fashion. 



The coefficients a(P) in (JT5j) are obtained by the following rule 2l(: For P — I, the 
identity, set a(I) = 1. For any other permutation decompose P into transpositions. For 
every transposition of the particles j and I write down a factor of —e %Q]l . For example, if 



P 



then one obtains the result 



where 





a(P) 



e iB it 



3 i(032+031+02l) 



c — — ki) 
c + i(kj — ki) 



or equivalently 



9ji = i log 



c + i(kj — ki 
c — i(kj — k{) 



-2 arctan 



(16) 



(17) 



(18) 



(19) 



where the branch of the logarithm or the arctan has not been specified yet. A more detailed 
description of how the coefficients a(P) are obtained can be found elsewhere j^ . 

The condition ()12j) determines only the form of the wave function. The allowed values 
for the set {kj} have to be determined by the periodicity conditions (JT3J) and (JHJ). In fact, 



these equations are equivalent to 

(-1) 



2l|: 



N 



N-l e ~ikjL 



exp(i^2e sj ), j = l,2,...,N, 



(20) 



where 8jj = 0. This set of coupled transcendental equations determines the wave vectors kj. 
However, in the present form it is still very cumbersome to work with, and it is customary 
to introduce new variables which exploit the symmetries of the problem. By solving ()20j) for 



kjL, one arrives at 



N 



kjL = 2Tcmj - 22 3 = !> 2 > 



N, 



(21) 



s=l 



for some integers {rrij}. The combination of ()21j) and (|19J) results in Lieb and Liniger's 
system of N — 1 coupled transcendental equations for the differences between the {kj}: 

rN c+i(ki-kj) 



(kj + i — kj)L = i Log 



c—i(ki—kj) 



n 



N 



c+i(k 



m=l c-i(k m -kj +1 ) 



+ 2t™„ j = l,...,N-l, 



(22) 
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where now the principal part of the logarithm is taken and the integers rij are discussed 
further below. Since the principal part of the logarithm is taken in (|22JI and not in (|2ip. rij 
is not necessarily equal to — rrij. Every state is uniquely defined by the set {kj}. The 
order of these is unimportant, since the particles are identical bosons. For real {kj} they 
can be ordered such that 

h<k 2 < ... < k N (23) 

which can be satisfied by choosing 

-2tt < $1(0 ji) < 0, (24) 
if j > I. The choice ()24)1 implies a one-to-one correspondence between the 9ji and the set 



{kj}, provided (j22j) have a unique solution for every set {rrij}, which can be justified |21|. 
When c varies continuously from — oo to +oo, the arguments of the logarithms in the 9ji also 
vary continuously. When the argument of any of the 9ji arrives at the discontinuity of the 
branch cut, the 9ji are continued analytically to compensate for the branch cut discontinuity. 
By taking the product of all N equations (|2~U|) one arrives at 

N N 



p = E^ = E 2 -T = 27r L- (25) 



rrin _ n 

271 

3=1 3=1 

p is the eigenvalue of the total momentum operator, which commutes with the Hamiltonian. 
Thus, the total momentum is quantized. Since the expression for p reduces to a sum over 
the rrij, it is an invariant of the 'motion' as c varies continuously from — oo to +oo. In 
particular, a state with zero angular momentum for c = will always have zero angular 
momentum, independent of the value of c. 

By inserting equation (j!5|) into ([lip, one finds that the full energy is given by the formula 

N 

Mc) = J2 k l ( 26 ) 

3=1 

It contains also the energy of the interaction between particles and can therefore not be 
considered ClS db purely kinetic energy. 

There is one more important thing to note about the possible {kj}. If the set {kj} is a 
solution, then the same is true for the set {kj}, defined by 

kj = kj + 2irn /L, (27) 



which shifts the total angular momentum by +27nioN/ L, where no is an integer. Therefore, 
only states in the central momentum strip have to be considered, namely the states with 

-nN/L <p< ttN/L. (28) 

To obtain a solution, the set of equations (}2*2*|) have to be solved together with and 
(J2HJ). Defining 

5, = (k J+1 - k,)L, j = l,...,N-l, (29) 

one changes variables from {kj} to {5j,p}. Technically this can be done by introducing the 
vectors 

k= (h,...,k N ) T , (30) 



S = (8i, ...,5 N -i,pL) 



(31) 



and the N by N matrix 



.4 



/-ll \ 
0-110 



0-110 
-1 1 

.. 1 1 1J 



V 1 1 ■■■ 

The transformation is then given by substituting 

1 



(32) 



(33) 



in (|22|). Even when N is large, the matrix A can easily be inverted. However, already for 
N = 4 the set of equations (|2*2*j) becomes so lengthy after this transformation that it makes 
no sense to give it here explicitly in terms of the new variables {Sj,p}. Still, it can always 
be obtained with any computer algebra program. In this work we used Mathematica [35( 
and in all numerical calculations we used the new variables {Sj,p}. 

The order (|23|) implies that Sj > when they are real. It has been shown that in the 
limit c — >• also c/8j — > [2]]. Hence, the arguments of all logarithms must approach unity 
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in this limit and thus are far away from the branch cut discontinuity. Therefore, if \c\ is 
small enough 

{ni, ...,n N -i} = {n$, ■■■,n° N _ 1 } (34) 

is certainly correct, where the n° denote the values of the rij for c = 0. The n° still remain 
to be specified. This can be done by considering the non-interacting ground state. 

D. General Consequences drawn from the non-interacting ground state 

In the absence of interaction the ground state wave function is simply a constant. This 
is equivalent to all kj and p equal to zero. Therefore, for c = one finds by using (JTSJl 

e*' 1 = -1, (35) 

and with the choice of the ranges of 6ji in (J2"2j) 

0;i = -tt, j>l, (36) 

and 

a{P) = 1 (37) 
for any permutation P. From (}2"T]) it follows that 

8j = 27r(m j+1 - rrij - 1) = 2vm°, j = 1, N — 1. (38) 

Hence all are zero for the ground state. In the present work this is the only state that we 
are interested in. The n° are unambiguously related to the rrij via (|38|) and (|25|). Therefore, 
they can be used equivalently to the rrij to classify all states. In this classification scheme 
the ground state is denoted by 

K,...,^} = {0,...,0} (39) 

Thus, the non-interacting ground state determines the values of the integers {n°, . . . , n^ r _ 1 }. 

The numerical computation for the interacting ground state can now be started with the 
values {ni, . . . , ri/v-i} = {0, . . . , 0} for the quantum numbers rij, provided |c| is sufficiently 
small (see previous subsection). However, the numerical effort can be reduced significantly 
by taking into account that all n° are equal for the ground state. For n° = const, it follows 
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from the symmetry of the set of equations (|22|) that if k is in {kj}, then so is — k 
This is equally true for repulsive and attractive interactions. Hence, for the ground state 
p is equal to zero. The order between the kj implies that k N+1 _j = —kj or equivalently 
<5jv_j = 5j for all j. Hence, the number of variables is reduced for even N to N/2 and for 
odd iV to (N — l)/2. This is also the number of equations that remain to be solved, as can 
be seen by substituting these equalities into (|2*2*j). and therefore we refer to this number as 
N eqs . Therefore, the problem consists now of solving N eqs coupled transcendental equations 
in N eqs unknowns. The numerical solution of these equations is discussed in section IIHI 

It is important to note that only the and not the rij may be taken for a unique 
classification of states since the latter may change when c is decreased or increased from 
zero onwards. In fact, for c > one finds numerically that the arguments of all logarithms 
in (|22|) move clockwise around zero on the unit circle when c increases. Hence, each of the 
rij has to be changed to rij + 1 when the argument of the j-th logarithm crosses the branch 
cut discontinuity of the principal part of the logarithm. For c < 0, it turns out that the 
ground state kj are purely imaginary and hence the arguments of all logarithms in f!22j) 
remain positively real. This implies that for attractive interaction the rij always remain at 
their values n° for zero interaction strength. Since the kj are purely imaginary for attractive 
interaction, it is customary to use the variables {otj,p} for c < 0, where 

aj = i6j, j = 1, . . .,N- 1. (40) 



E. Perturbation theory 

For sufficiently weak repulsive or attractive interaction perturbation theory should also 
be applicable. Since the non-interacting ground state is just a constant, it is easy to derive 
the energy to first order. By treating the whole interaction potential 

N 



2c^2d(xi- Xj) (41) 



Kj 

as a small perturbation of the unperturbed ground state with the wave function 

N 



if)(xi,...,x N )=[-=) , (42) 
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one finds by using non-degenerate first order perturbation theory that the ground state 
energy per particle is given by 

4 1} (c) c(N-l) 



N 



(43) 



III. NUMERICAL SOLUTION 



In this section we discuss the numerical solution of the set of coupled equations 
starting from the non-interacting ground state. By solving (J22j) . it is assumed that the 
ground state takes on the form of a Bethe-ansatz wave function. For repulsive interaction 



the validity of the Bethe-ansatz has been proven rigorously 
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3- 



However, for attractive 



interaction there is no such proof to our knowledge. As discussed in section |Hj for two and 



three bosons there are known cases 
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3ll | in which the form of the wave function is not of 
the Bethe-ansatz type for certain critical, attractive interaction strengths. However, for the 
ground state this assumption is justified further below by using an argument similar to that 
of Muga and Snider |3l| . For all computations we used the Mathematica function FindRoot 
35| which makes use of the Newton-Raphson algorithm. All calculations were carried out 
on a standard PC. 

We define the residual error A res £ls the sum over the absolute values of the differences 
between the left- and the right-hand sides of each of the equations (|2"2"|) 

Neqs 

A res = | A,- 1, (44) 

3=1 

where Aj is the difference in the j-th equation. Our goal is a residual error of 

A res < 1(T 9 (45) 

in all computations. Thus, we ensure that the wave vectors kj are accurate up to the eighth 
digit. We found that for repulsive interaction this goal is fairly easy to achieve, whereas 
the attractive case proved to be a lot more problematic. In the numerical calculation, we 
started from the non-interacting ground state, and increased/decreased c stepwise. For very 
weak interaction and small particle numbers the FindRoot function is not too sensitive to 
the initial guess. As an initial guess for the following sets of {5j} or {aj} we used the three 
previously calculated sets and extrapolated to the next one. 
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The computed ground state energies per particle are depicted in Fig. [T] as a function of 
the interaction strength c. It is seen that the energies of the repulsive ground state solutions 
are smoothly connected at c = to the energies per particle of the corresponding Bethe- 
ansatz solutions for attractive interaction. The smooth connection of the energies is by no 
means trivial, as we shall show now. As mentioned in section |HJ the Sj are purely real 
for repulsive interaction, but purely imaginary for attractive, motivating the redefinition of 
variables (}4T)j) . This can be seen for the case of fifteen particles in Fig. |21 which we shall 
discuss representatively for all other particle numbers. Since the 5j are real for c > and 
imaginary for c < they can not be smoothly connected at c = 0. It is therefore quite 
surprising that the energy, which is a function of these variables, is smoothly connected at 
c = 0. 

For repulsive and attractive interactions it can be seen that all 5j and all otj respectively 
converge to zero when |c| — > 0. Together with p = this implies that the non-interacting 
ground state for which all kj are zero is approached from either side. 

In the repulsive case all 5j start from zero and begin to spread when c is increased from 
zero onwards. After a maximal spread the Sj start to degenerate successively and saturate 
to their value at infinity, namely 27r jsij ] . 

For attractive interaction the variables ay are used. The atj behave similar to the repulsive 
case in the sense that they first start from zero, spread and then degenerate, but there is no 
saturation for strong attractive interaction. These findings are equally true for all particle 
numbers studied in this work. 

The smooth connection of the energies and the continuous connection of the kj at c = 
prove the validity of the Bethe-ansatz numerically for the ground state of up to fifty bosons. 
However, we suspect that a complete set of states can be obtained for any particle number 
from a Bethe-ansatz wave function in the sense explained in section ITT1 

Numerically, the situation is very different for repulsive interaction and attractive inter- 
action. The degeneracy sets in much earlier for attractive than for repulsive interaction. 
This is important for the following reason. We found that the numerical calculation tends 
to break down, if the relations 

<5i >5 2 > ... > 6 Neqs , (46) 
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cti > a 2 > ... > a Neqs (47) 

for repulsive and attractive interaction, respectively, are not fulfilled at all times. This 
illustrates the importance of a proper initial guess. The graphs of the {5j} and {a,} for 
other particle numbers are very similar to Fig. 121 We discuss the effects that arise due to 
changing N in section HVl and section IVl 

Due to the degeneracy of the 8j and the otj we were forced to perform all calculations 



with high precision numbers. The Mathematica option WorkingPrecision |35J] allows for 
computations with numbers of arbitrary precision. For repulsive interaction we found that it 
is sufficient to work with a number precision of no more than 10 -20 for all variables involved 
in the computation, at least for up to fifty particles and c < 1500. Then, the smallest 
difference between the 8j, £24 — #25, is still greater than 10 -6 . For attractive interaction on 
the other hand we found that we had to use numbers of extremely high precision in the 
calculations for fifty particles, even for |c| < 1. Then, the difference «24 — 0^25 is as little as 
10~ 85 . For fifty particles we could achieve our goal of a residual error of less than 10~ 9 only 
by using numbers with 90 digits! However, the fact that all computations were carried out 
on a standard PC proves that the possibilities of our approach are still far from exhausted. 



IV. RESULTS ON REPULSIVE BOSONS 

A. Exact finite N solution for repulsive interaction 

By solving the set of coupled equations for c > 0, as described in section 11111 we 
find for each c the corresponding set {5j,p}. Using (J33|) . one obtains the corresponding set 
{kj}. In Fig. El the set {kj} is depicted for the case of fifteen bosons. Starting from the non- 
interacting ground state, where the wave function is a constant and all kj are zero, the kj 
spread for c > 0. In the limit of very strong interaction we recover Lieb and Liniger's result 
that all kj become constant and equally spaced with a separation of 2tt between adjacent 
kj. The expression for the energy (ffilj) then implies that the energy per particle saturates 
to a finite value when c —>■ 00. For other particle numbers the situation is much the same, 
only that the larger the number of particles, the later this saturation sets in. 



15 



B. The relation to the Tonks-Girardeau limit 



The saturation for strong interaction can be quantified by considering the energy dif- 
ference to the saturation energy From a physical argument it is clear that the saturation 
energy has to coincide with the energy per particle of the TG gas which is given by 2^]: 

~ 3L* ' (48) 

where L is the length of the ring. In Fig. EJthe exact energies per particle on a ring of unit 
length are plotted as a function of c and the TG energy is indicated on the right border of 
the graph for different particle numbers. It can be seen clearly that the larger the number 
of particles is, the stronger the interaction has to be to reach a given fraction of the TG 
energy. 

We would like to answer the following question quantitatively: What is the smallest c for 
which the TG limit is less than a certain percentage r away from the exact energy of the 
system? This relative deviation is given by 

and thus we are looking for the c values c r , for which 

' AE 



E 



(C.) = r. (50) 

L,TG 



This relative deviation may also be considered as the error in the energy introduced by 
approximating the N bosons on a ring at finite c by the TG limit. 

First, we discuss the scaling properties of {^f) LTG ( c ) with respect to L. By using (JJJ) 



and (SHI) one finds that 



f AE\ , . / AE\ . L' 

c ) = hr (c-r). (51) 



V E J L i TG \ E J L TG L 

Relation (jSTjl allows to calculate the values c r for a ring of length L', if c r is known for a 
ring of length L. The result is shown in Fig. |3J For illustration purposes we plot the c r 
values for a constant density p = N/L = 1, rather than a constant length of the ring. To 
get the corresponding c r values for a ring of length L = 1, one simply needs to multiply 
each c r value by N. Fig. El shows that the c r are almost independent of iV for any given r, 
when the density is constant. For a ring of length L, this implies that the corresponding c r 
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depend almost linearly on the number of particles. The slopes of the curves for the c r on a 
ring of unit length as a function of N are exactly the c r in Fig. |S| The inset shows how little 
the deviation from a constant slope is, even for r as large as 95%. When r is decreased this 
deviation becomes even smaller. The fact that the TG limit is reached for larger c-values 
with increasing N can therefore solely be attributed to the increase in the density by adding 
more particles. Using ()48|) and Fig. El one can even obtain a rough estimate to the exact 
ground state energy for any particle number, simply by extrapolating for the desired N on 
the graph to c, multiplying c times N and solving (|50|) for E L= i(c). It is surprising, how 
simple the relation between the system for finite c and the TG limit is at constant density. 
There is almost no dependence of the c r on the particle number, especially when r is small 
and N is more than just a few particles. 

C. The relation to the thermodynamic limit 

Now we address the question, how far away the thermodynamic limit [21] is from the 
finite N solution. In the thermodynamic limit N, L — > 00 with p = N/ L = const, and 
the energy of this solution differs from the exact energy of the finite N system, even if the 
densities are the same. The energies E p= i^dl/N of this solution can be found elsewhere 
(3^). For large particle numbers, N ^> 1, it is expected that the thermodynamic limit energy 
gives a good approximation to that of the finite N system, provided the densities are the 
same. For small N, however, the finite number of particles should play an important role. 
We would like to answer the following questions. Firstly, how large does N have to be in 
order that the thermodynamic limit solution provides a good approximation to the finite 
system, and secondly, what is the quality of this thermodynamic limit approximation, when 
c varies for a given N? 

To answer these questions quantitatively we consider the energy of the thermodynamic 
limit solution for the density p = N/L = N and compare with our results on a ring of 
unit length. This ensures that we compare the two systems at the same density. We define 
the relative deviation of the thermodynamic limit energy from the energy of the finite N 
solution 



This quantity can also be considered as the error which is introduced by using the energy 




(52) 
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of the thermodynamic limit solution instead of that of the finite N solution. Again, we first 
discuss the scaling properties of this relative deviation when L varies. The ground state 



energy in the thermodynamic limit can be written as 



2l|: 



E p ,tdl(c) = Np 2 e{-\ (53) 



where e(x) is a monotonically increasing function which is tabulated in [36]. Therefore, we 
find the same scaling behaviour as for the TG limit (|51|). as follows by using ((7|) 

if) <^(f) «£• <«> 

\ ^ / L',TDL V ^ J L.TDL ^ 

The relative deviation LTDL (c) is shown in Fig. El for = 1. At first sight, the result 
is quite surprising. For all particle numbers the relative deviation is largest for c = 0, but 
decreases rapidly before saturating to a finite value for infinitely strong interaction. However, 
one can explain this behaviour by considering the limiting cases c — > oo and c —>■ 0. 

We begin with the limit c — > oo. In the thermodynamic limit the impenetrable boson gas 
energy per particle is given by [22]: 

E TV 2 

= yP 2 , (55) 

whereas (14*%)) is the c — > oo limit of the finite N system. Substituting L = 1 in PH|) and 
p = N in (JST)|) results in 

^L.™ ^ (56) 

which is exactly the tendency that can be seen in Fig. IHJA 

For zero interaction strength the behaviour can be explained by treating the full interac- 
tion between particles as a small perturbation to the non-interacting ground state for finite 
N, which results in (}4*3*|) with L = 1. The corresponding 
particle in the thermodynamic limit for small c is given by 

E p= n,tdl{c) 



eading term of the energy per 
3: 



N 

Hence, the relative deviation {^r) L - lTDL for c — > becomes 

fAE\ c ^o 1 



cN. (57) 



\ ; / V i (58) 

\ e j L=hTDL Jy - i 



as can be seen in Fig. U3. 
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An interesting consequence of (|54|) is that for given N and L, the relative deviation 
(^F)z tdl ( c ) * s obtained by the value of the curve for the same N in Fig. El evaluated 
at c = Lc. Varying the length of the ring only changes the point of evaluation in Fig. El 
Consequently, neither (|5fi|) nor (|58j) depend on the length of the ring. Obviously, for small 
particle numbers the thermodynamic limit approximation is never a good approximation, 
no matter how strong the interaction or how large the size of the ring is. 



D. Quality of the thermodynamic limit approximation 

In the context above, the following question arises naturally. Given a certain interaction 
strength c and a ring of a fixed size L, how many particles are at least necessary for the 
thermodynamic limit approximation to be accurate to a certain percentage, for instance 
1%? We denote this number by N\%. The result is shown Fig. For weak interaction 
strength the number of particles has to be much larger than for strong interaction. Although 
this is counterintuitive, one should remember that weak interaction is the regime where the 
thermodynamic limit approximation is worst, see Fig. El Equivalently, for a given particle 
number, the size of the ring has to be much larger for weak interaction strength than for 
strong interaction, if the thermodynamic limit approximation is to be used for the description 
of the experiment. It is no coincidence that the curves in Fig. 0all look the same. From 
(JH4"J) it follows that if one changes L — > V and simultaneously c — * p-c, then the deviation 
(^r) LTDL remains unchanged. Therefore, the curve, e.g., for c = 0.1 can be obtained from 
the curve for c = 1 by multiplying all L values times ten. This permits to obtain curves for 
all values of c. 

From (Jo7)|) it follows that for iV < 10, the relative deviation {^) L - 1TDL does not drop 
below one percent for any c. Moreover, from (JH4"j) it follows that this can not be compensated 
for by making the ring larger. We conclude that the thermodynamic limit approximation 
can never be accurate to 1%, if the number of particles is less than 11. Similarly, it follows 
from that only for iV > 101 the relative deviation (^P) LrDL is always less than 1%. 
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V. RESULTS ON ATTRACTIVE BOSONS 



A. Exact finite N solution for attractive interaction 



In this section we discuss the results obtained by solving the coupled equations ({22} for 
c < 0, as described in section IIHI Again, the kj are shown for fifteen particles in Fig. |S1 
Since the 5j are purely imaginary, only the imaginary part of the kj is plotted. Similar to the 
repulsive case, the kj start to spread when |c| is increased. However, for attractive interaction 
there is no saturation when |c| — > oo and the kj keep on spreading. As can be seen from Fig. 
121 the aj and therefore also the kj grow virtually linear with |c|, already for comparatively 
weak attractive coupling, implying a quadratic dependence of the energy on c. This is the 
region in which the aj are practically degenerate. The larger the number of particles, the 
earlier the degeneracy of the aj sets in. Fig. EJA. shows the energy per particle for different 
particle numbers as a function of |c| on a ring of unit length. A linear dependence in the 
vicinity of c = is followed by a quadratic decrease for stronger interaction. Consistent with 
the perturbation theory results (|43|) . the energy per particle for larger particle numbers is 
below that for smaller particle numbers. 



B. The relation to the system on an infinite line and to mean-field theory 

While in the repulsive case it is possible to compare the energy to either the TG limit or 
to the thermodynamic limit, neither of these exist in the attractive case. This follows from 
the behaviour of the exact solution, as discussed below. For attractive interaction we found 
that a different system is closely related to the N bosons on a finite ring. Namely, this is 
the system of N bosons subject to attractive 5-function interaction on an infinite line. This 
system has been solved exactly, see and references therein. The energy per particle is 
given by 

= ~<?{& - !)• (59) 

Although no periodic boundary conditions are imposed in this case, it is possible to think of 
this system as the limit L — > oo of the system on a ring of length L. However, the relation 
between these two systems is far less evident than the relation of the repulsive system to the 
impenetrable boson gas and the thermodynamic limit solution. Therefore, we will motivate 
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this relation in a simple mean-field picture before we begin with the comparison. We briefly 
review the main mean-field results for the attractive case. 

In the case of Bose-Einstein condensates the standard mean-field approximation is the 
so called Gross-Pitaevskii (GP) equation. The GP equation has been solved ana" 
the problem that we consider, namely the one-dimensional condensate on a ring 
Since the interaction strength appears in the GP equation only in combination with the 
particle number, the interaction strength can be parameterized by the new parameter 

Thus, increasing the particle number or the interaction strength are equivalent in this ap- 
proximation. The ground state wave function is simply a constant as long as \G\ is less 
or equal to a critical |G C r,L|, where G ct> l = — f for a ring of length L. It is an angular 
momentum eigenfunction with zero angular momentum. As soon as GjG cr ^ > 1 a second 
solution appears which is lower in energy than the constant solution, but is not an angular 
momentum eigenfunction, since it localizes at some arbitrary point on the ring. For very 
strong attractive interaction it is essentially zero everywhere, except for the position around 
which it started to localize, much like a 5-function. In this mean-field picture it is clear that 
the length of the ring loses its importance in the strong interaction limit, since the wave 
function is localized on a small fraction of the ring and hence does not "see" the finite size 
of the ring. It should therefore make no difference how large the ring is. Furthermore, this 
effect should become important when G/G cr ^ > 1. Of course, this mean-field picture is 
oversimplifying, since the true many-body wave function must be an eigenfunction of the 
angular momentum operator, but it captures some aspects of the true situation, as we shall 
show below. 

We return to the discussion of the results. Fig. EfB shows the energies per particle for 
different particle numbers as a function of the parameter G/G cr ^ for L = 1. Since G contains 
a factor of N — 1, for any constant G the energies per particle for larger N are now above 
that for smaller N. Also the GP ground state energy is depicted. 

Now, we consider the difference between the exact energies per particle on an infinite 
line and on a finite ring (El =00 (c) — El(c))/N. If the mean-field picture discussed above is 
correct, a sharp decrease in this energy difference is to be expected when the GP solution 
starts to localize. This should take place at G — G cr ^ and therefore the corresponding 
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critical c values c cr = should be proportional to 1/(JV — 1). Fig. ITUK shows that this 

is indeed the case. For all particle numbers the energy difference takes on a maximum, just 
before G = G cr ^ and decreases rapidly afterwards. 

We now consider the relative and not the absolute deviation of the energy of N attractive 
bosons on an infinite line from the energy of those on a ring of length L. This deviation is 
given by 

Similar to the limits to which we compared in the repulsive case, this is the error which 
is introduced by approximating the ring by a line of infinite length. First, we discuss its 
scaling behaviour when L varies and then its properties when considered as a function of 
G/G crt L rather than c. Using (J7|) and the quadratic dependence of on c one finds that 
the relative deviation on a ring of length L is related to the relative deviation on a ring of 
length L' via the equation 

The relative deviation J (c) for attractive interaction scales exactly in the same way 
under a change of L as the relative deviations {^jjr) LTDL ( c ) an d (xliTG ( c ) ^ or re P u l s i ve 
interaction, although the limits considered are completely different. It follows from ()62j) that 



the relative deviation {^^j (c) for a ring of length L is simply the relative deviation on 
a ring of unit length, evaluated at c = Lc. 

How are the relative deviations ('^f~ L ) ( c ) an d ( c ') on ^ wo r i n g s °f length 

L and V related when c and d are chosen to satisfy G/G^l = G' /G cr ^l In this case 
the interaction strengths on both rings equal the same fraction of the critical mean-field 
interaction strength on each ring. The relation G/G cr> L = G' jG cr ^ is equivalent to d = Cjj 
which implies (c) = \-^r-) ( C Q when substituted into The relative deviation 

V L J OO \ L J OO 

on rings of different lengths is exactly the same, if the interaction strengths are the same 
fraction of the critical mean-field interaction strength on each ring. As a function of G/G cr) L, 
it is therefore independent of L. This relative deviation is depicted in Fig. HUB. Its absolute 
value decreases linearly for GjG cr ^ < 1 and exponentially for GjG cr ^ > 1- It can be seen 
clearly that for interaction strengths that are about twice as strong as the critical mean- 
field interaction strength the energy of the system of iV bosons on a ring is given to a good 
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approximation by the energy of the system on an infinite line. This is equally true for all 
particle numbers studied. Moreover, the factor 1/L in G CT ^ implies that the larger the ring 
is, the smaller c has to be for (|59j) to be a good approximation to the exact ground state 
energy of the finite ring. 

The fact that (E L=00 (c) — E L (c))/N and the relative deviation {^^j (c) converge to 
zero proves that the simple mean-field picture of a localizing wave function reproduces at 
least some of the true physical behaviour. Interestingly, (E L=OQ (c) — E L )(c)/N is greater 
than zero, implying that the ground state of N attractive bosons on an infinitely long line 
lies above the ground state energy of those on a ring of finite size. Clearly, this is against 
the physical intuition. Energy levels are normally lowered when the distance between the 
confining boundaries of a system is increased. However, since periodic boundary conditions 
are used, there are no confining walls and this anomalous behaviour must be attributed to 
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3 11 ] . This anomalous 



the periodic boundary conditions, as was already pointed out earlier 
behaviour is obviously present for any particle number. 

The weak-interaction behaviour of (^ij^J ( c ) can also be found by using (|43 p .([59 |) and 
(E2J): 

(AE L \ . . c -»o 1 | I r / at , -I \ i ^ G N+l 
{^r)J C) -T2^ L ^ N+1) - 1 = -2G^N^T- 1 - (63) 

It can be seen in Fig. 110b that the perturbative expression (J53)l gives a good approximation 
to (c) for large particle numbers in the whole region G /G cr ^ < 1, whereas for small 

particle numbers the validity of the approximation is very limited. This can be understood 
by remembering that the c values corresponding to G cr ^ are proportional to N — 1. For large 
N the critical c are sufficiently close to zero and perturbation theory becomes applicable. 
Therefore, only for large N the expression gives a good approximation to (^^^j (c) 
in the whole region GjG CT ^ < 1- 

To conclude the discussion in the mean-field picture, it has to be mentioned that on an 
infinite line not only the exact ground state solution, but also the GP solution is known 
analytically. Their energies are related by 

C2=^-^, (64) 

implying that for N ^> 1 and/or G/G cr> L ^> 1 even the GP energy E^J2> gi yes a good 
approximation to the exact energy of the finite system. This is in contrast to the repulsive 
case, where the energy of the corresponding GP solution diverges for large G and therefore 
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does not reproduce the true physical behaviour, i.e. the TG limit [11]. As mentioned earlier, 
the GP approximation unfortunately breaks the rotational symmetry of the Hamiltonian. 
However, the situation can be remedied. Restoring the broken symmetry provides a ground 
state wave function of lower energy than the GP solution and unique many-body properties, 
see 

Since the energy of the system of iV bosons on an infinite line ()59|) provides a good 
approximation to the energy per particle of the system on a finite ring, provided |c| is not 
too small, it allows to investigate the possibility for different limits in the attractive case. 
The energy per particle diverges with c 2 for strong interaction and therefore there is no 
attractive equivalent to the TG limit. Similarly, the thermodynamic limit does not exist 
since the energy per particle is always below (J59|) which is independent of the size of the 
ring and diverges with the number of particles as iV 2 — 1. 



C. The importance of the finite length of the ring 

As we have shown above, the system of N attractive bosons on a ring is closely related to 
the system of N attractive bosons on an infinite line. We found that the energies of these two 
systems are the closer, the larger the number of particles and the stronger the interaction is. 
Now, we would like to address the following question: Given a ring of length L and a fixed 
interaction strength c, how many particles are at least necessary that the energies of the 
two systems mentioned before differ by no more than one percent. We denote this particle 
number by Ni%, keeping in mind that Ni% has a different meaning for repulsive interactions. 
Thus, we are looking for those particle numbers for which the relative deviation j (c) 
defined in ()61)1 is one percent. It is possible to think of N\% as the particle number from 
which onwards the finite ring may be approximated by an infinite line within error bounds 
of one percent. The result is shown in Fig. EI For all particle numbers (c) drops 

below one percent, if either L or c is large enough. This implies that it is possible for 
any particle number to approximate the finite ring by an infinite line either by making the 
interaction stronger or the ring larger. From (p)2|) it follows that if one changes L — > V and 
simultaneously c — > jjc, then the deviation (c) remains unchanged. The curves can 

therefore be shifted horizontally in the sense that a change of c — > \c is accompanied by a 
change L — > xL. 
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The fact that the impact of the finiteness of the ring vanishes for strong interaction 
manifests itself also in the degeneracy of the aj which are depicted in Fig. By considering 
the differences between the <x,- as an effect which is due to the finite size of the ring, it is 
possible to derive an asymptotic relation between the energy of the system on an infinite line 
and that on a finite ring. For any fixed value of c it is possible to choose a length of the ring 
Lq such that G/G cr ^ 3> 1 for any L > L . The acj are then virtually degenerate although 
they satisfy (JTTj) at all times. This can be seen in Fig. 03, where G/G C r,L = 1 corresponds 
to |c| ~ 0.7. To analyze this degeneracy we set <x,- = a = const for j = 1, . . . , N — 1 All 
wave vectors kj are then equally spaced. We calculate the energy of the system in this case, 
keeping in mind that for the ground state always —A; is in the set {kj} when k is. By using 
(I29j) and (J4(J|) one finds that the energy (J26|) for an odd number of particles is given by 

■Y-l 

N — 2 — 

e =E*? = 2 £h^> 2 = ~h (iT n(n2 ~ l) - (65) 

This is exactly the ground state energy of N bosons on an infinite line ()59jl . if we put 

a = \c\L. (66) 

An analogous calculation for even N leads to the same expression. We have thus found 
the asymptotic c-dependence of the ground-state atj for c — ► — oo and any particle number. 
The resul t ()66|) for iV particles is consistent with the result of Muga and Snider for three 
particles 3l|. In principle, it is also possible to calculate a correction to this asymptotic 
c-dependence, for instance, by applying the techniques developed by Muga and Snider (31 1. 
but we refrain from doing so due to the complexity of the equations f!22[) for large particle 
numbers. 



VI. CONCLUSIONS 



In this work we have presented the exact ground state solutions of the many-body 
Schrodinger equation for up to fifty bosons on a ring subject to pairwise 5-function interac- 
tion. By employing the Bethe-ansatz for attractive interactions we have proven numerically 
that the Bethe-ansatz provides ground state solutions not only for repulsive interaction, 



where this is well known to be the case 
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33j , but also for attractive and all particle num- 
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bers studied. We suspect that a complete set of states can be derived from a Bethe-ansatz 
for any particle number and any interaction strength. 

Our results show that the repulsive and the attractive case exhibit fundamentally different 
behaviour. While in the repulsive case all wave vectors are purely real, they are purely 
imaginary in the attractive case, implying a jump in the derivative at the point of zero 
interaction strength. For very strong repulsive interaction the wave vectors converge to 
finite values which are equally spaced, in agreement with the result of Lieb and Liniger 



2l|. For strong attractive interaction on the other hand the wave vectors do not converge 



to finite values. We have derived an asymptotic relation for the wave vectors in this limit. 
Asymptotically, the wave vectors for attractive interaction are then equally spaced, while 
their absolute distance keeps on growing with increasing interaction strength. 

For repulsive interaction our exact data are compared in some detail with the energies 
determined in the thermodynamic limit and in the Tonks-Girardeau limit. For attractive 
interaction, where these limits do not exist, we have compared the exact energies with those 
analytically known for a system on an infinite line. 

In the repulsive case we have shown that the energy of Lieb and Liniger's well known 
thermodynamic limit solution can differ substantially from that of our exact solution, espe- 
cially when the number of particles is small and the interaction is weak. In detail, we found 
that the approximation of the finite system by the thermodynamic limit solution never re- 
produces the exact energy within an accuracy of one percent, if the number of particles is 
less than eleven. Moreover, our investigation has revealed that this thermodynamic limit 
approximation is accurate up to one percent for all repulsive interaction strengths only when 
the particle number is at least as large as 101. The results obtained allow one to conclude 
when the thermodynamic limit approximation is a good approximation and when finite N 
effects have to be taken into account. This is not only of academic, but also of practical 
importance since in recent experiments on Bose-Einstein condensates only a few dozens of 
atoms could be studied, see Jfjj and references therein. Hopefully, our results will be helpful 
to improve the description of these experiments. 

We have also investigated the system for finite c in the light of the Tonks-Girardeau 
limit. For strong repulsive interaction the energies of our exact solutions approach the 
Tonks-Girardeau energies. Our analysis proves that very strong repulsive interactions are 
necessary to approximate the exact energy by the Tonks-Girardeau energy within some well- 
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defined error bounds. However, we found that the convergence towards the Tonks- Girardeau 
limit is virtually independent of the number of particles when the particle density is held 
constant. The relation revealed between the Tonks-Girardeau limit and our exact solution 
allows to estimate the exact ground state energy of the system for any particle number and 
any repulsive interaction strength. 

In the attractive case we have related the system on a ring of finite size to that on an 
infinite line. The ground state energies of these two systems are found to be essentially 
identical when the interaction is much stronger than a certain critical interaction strength. 
This critical interaction strength originates from mean-field theory and we could explain the 
strong relation between these two systems in a mean-field picture. 

Interestingly, the energy of the system on the finite ring is always below the corresponding 
energy of the system on an infinite line. This anomalous behaviour was already found earlier 

21L |3l| and was found to be present for any other particle number 



for two and three particles 
studied in this work. 

For all interactions we have given explicit bounds on the minimal number of particles and 
the minimal size of the ring for the reference systems discussed to be good approximations 
to the exact solution. 

Although the ground state for repulsive and attractive interaction is obtained by solving 
the same set of coupled transcendental equations the numerical effort for solving these 
equations differs substantially in the two cases. While the system remains reasonably stable 
for repulsive interaction, the attractive case requires the use of very high precision numbers 
already for weak interactions. 

Finally, we stress that the excitation spectrum of the system which is of high interest by 
itself can also be obtained with the approach presented in this paper. 

The explicit numbers of the exact ground state energies of the finite N system studied 
in this work can be found on the Internet j^J]. 
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FIG. 1: (color online) Exact ground state energies per particle as a function of the interaction 
strength c. The Bethe-ansatz was used to compute the solutions to the problem of N bosons on 
a ring of length L for repulsive c > and attractive c < interaction strengths. It can be seen 
that the ground state energies for repulsive interaction are smoothly connected at c = to those of 
the Bethe-ansatz solutions for attractive interaction with the same particle and quantum numbers 
(see text). This proves the validity of the Bethe-ansatz also for attractive interaction for all cases 
studied. The ring was taken to be of unit length L = 1. 
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FIG. 2: Lieb and Liniger's set of coupled transcendental equations is solved by the set of variables 
{5j}. For the repulsive ground state all Sj are real, whereas they are purely imaginary for the 
attractive one. This motivates the use of the variables ctj = iSj for attractive interaction. Depicted 
are the sets {5j} and {aj} for fifteen bosons on a ring of unit length. In the absence of interactions 
all Sj are zero. A: Repulsive interaction. Starting from zero the Sj begin to spread and are 
well separated from one another before approaching their limiting value of 2n. B: Attractive 
interaction. The aj spread, but start to degenerate already for comparatively weak interaction. 
This complicates the numerical solution significantly. 
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FIG. 3: Wave vectors kj for fifteen repulsive bosons on a ring of unit length. Only for very strong 
interaction the wave vectors are practically equally spaced with a separation of 2ir between adjacent 
kj, also see Fig. 
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FIG. 4: (color online) Ground state energies per particle for different N on a ring of constant 
length (L = 1) for repulsive interaction (c > 0). For c — > oo the energy per particle converges to 
that of the Tonks- Girardeau (TG) limit. On the right border of each graph the energy of the TG 
limit is indicated for each particle number. For larger particle numbers the TG limit is approached 
for larger values of c. 
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FIG. 5: Measure for the TG energy to approximate the exact energy. Given N bosons on a ring of 
variable length, at the interaction strength c = c r the relative difference of the TG energy to the 
exact energy equals the fraction r (see definition in the text and equation (|49jl). For this graph the 
density p = N/L was held constant, p = 1. Surprisingly, for a given r the values c r are practically 
independent of the number of particles when p is a constant. Even more surprising is the fact that 
this is equally true for small and large values of r. Top to bottom: r = 5%, 10%, 25%, 50%, 95%. 
The inset shows how little the deviation from a constant is, even for r as large as 95%. The 
corresponding c r for a ring of length L = 1 can be obtained by multiplying each c r value by N. 
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FIG. 6: (color online) Measure for the thermodynamic limit to approximate the exact solution. 
Lieb and Liniger's solution in the thermodynamic limit is often used as an approximation to the 
true finite N solution. This introduces an error which depends on the number of particles, the 
interaction strength and the length of the ring. In the limits c — > and c — > oo the error in the 
energy can be calculated by using perturbation theory (see text) and the TG expression for the 
energy (see text). The formulas obtained for infinite and zero interaction strength, see Eqs. (|56|) 
and (|58j) respectively, imply that this error is never less than 1% if N < 11. Similarly, only for 
N > 101 it is always less than 1%. These results do not depend on the length of the ring, (see 
text). Curves from top to bottom: N = 2, 3, 4, 5, 11, 15, 25, 50. Length of the ring: L = 1. A: The 
error introduced by using the thermodynamic limit solution instead of the solution for finite N is a 
monotonously decreasing function of the interaction strength and the number of particles. In the 
limit c — > oo the error approaches the value N -l_ 1 ■ B: For c — > the error converges to the finite 
value jr^r in agreement with perturbation theory. 
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FIG. 7: Depicted are those N for which the error in the energy introduced by using the thermody- 
namic limit solution instead of the finite iV solution is less than one percent. These N are denoted 
by Ni%. Equation (|54[) implies that the curves can be shifted horizontally, in the sense that a 
simultaneous change of L — > xL and c — > does not change Ni%. It can be seen that there is 
no Ni% < 11 since for N < 11 the error of the thermodynamic limit approximation never drops 
below one percent (see also Fig. |BJ). 
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FIG. 8: The ground state wave vectors kj of fifteen bosons for attractive interaction on a ring 
of unit length. All kj are purely imaginary. In contrast to the repulsive case, see Fig. |31 there 
is no saturation for strong attractive interaction. For strong attractive interaction all kj grow 
(approximately) linearly and are (almost) equally spaced, already for comparatively weak interac- 
tion. The interaction strength at which this transition to the (almost) linear dependence occurs is 
proportional to 1/(N — 1) (see text). In this graph the critical mean-field interaction strength is 
c « 0.7. 
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FIG. 9: (color online) Energies per particle for attractive interaction. The system resembles the 
system of N attractive bosons on an infinite line. This can be explained in a mean-field picture 
(see text). A: Energies per particle as a function of |c|. With increasing particle numbers the 
energy per particle drops dramatically. B: Energies per particle as a function of G/G cr ^L, where 
G = 2c(N — l)/(27r) and G cr ,L = —tt/L is the critical and fixed mean-field interaction strength. 
The GP energy is also shown (dashed line). For large N the GP energy approaches the exact 
energy. Length of the ring: L = 1 
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FIG. 10: (color online) On the ability of N attractive bosons on an infinite line to approximate 
the system of N bosons on a ring. A: Shown is the difference of the exact energies per particle 
{El=oo — El) / N of N attractive bosons on an infinite line and on a finite ring [L = 1) as a function 
of G/G cr ,L- This energy difference is counter intuitively positive, implying that the energy on the 
ring of finite size is below that of the system on an infinite line. It starts to decay to zero for all 
particle numbers at approximately the same value of G/G cr> i,. B: Shown is the relative error which 
is introduced by using the energy of the system of N attractive bosons on an infinite line L = oo 
instead of the energy of the system on a ring of finite size. As a function of G/G cr ^i this graph is 
the same for all values of L. This relative error decreases linearly for weak interaction, followed by 
an exponentially decaying tail. 
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FIG. 11: Depicted are the values of N - denoted by N ± % - for which the error introduced by using 
the energy of N attractive bosons on an infinite line instead of the energy of N bosons on a ring of 
size L is less than one percent. This error drops below one percent for all particle numbers once L 
is large enough. The curves can be shifted horizontally in the sense that a simultaneous change of 
L — > xL and c — > does not change iV^. Note that the meaning of is different for repulsive 
and attractive interaction, since the respective exact results are compared with different limiting 
situations. 
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